Purely electronic THz polarization in dimer Mott insulators 
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We theoretically discover purely electronic polarization modes in THz frequency region in dimer 
Mott insulators k-(BEDT-TTF)2X. The unusual low-frequency modes arise from the coupling be- 
tween the oscillation of intradimer electric dipole moments and that of alternating interdimer bond 
orders. These collective motions play an important role in the dynamical dielectric properties of 
the dimer Mott insulators. Near the phase boundary of the dimer Mott transition, the ferroelectric 
ground state is realized by introducing electron-lattice coupling. 
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Low-dimensional strongly correlated electron systems 
with a quarter-filled band have recently attracted much 
attention. BEDT-TTF salts are typical two-dimensional 
systems with anisotropic triangular lattices, and the 
strong correlation and frustration act cooperatively to 
generate exotic phases. Among the BEDT-TTF salts, re- 
(BEDT-TTF) 2 X (X: a counter anion) exhibit the Mott 
insulator phase [3-0, exotic superconducting phase [f|, 
spin liquid phase \ffi%, and so on. The BEDT-TTF 
molecular lattice is distorted to form dimers in k-(BEDT- 
TTF) 2 X. Hybridized two sites can be effectively treated 
as a single site, and the valence band is regarded as a 
half-filled one. As a result, k-(BEDT-TTF) 2 X exhibits 
the Mott insulator phase if the effective on-site Coulomb 
energy is large enough. The state is called the dimer 
Mott insulator. 

Recently quite interesting charge excitation phenom- 
ena have been observed in the dimer Mott insulators: 
photoinduced phase transition to a metallic state [Toj j . 
and anomaly in dielectric properties [ll| . In particular, 
the latter suggests that low-energy charge excitation ex- 
ists in the dimer Mott insulators in contrast to the case 
of the Mott insulators. 

In this letter, we show for the first time that the dimer 
Mott insulator has THz electric polarization, which 
purely arises from the electron dynamics. The ferroelec- 
tricity which arises from the electron dynamics has been 
observed in some materials, and these ferroelectrics have 
unconventional origins such as charge order generation 
We propose another new origin of ferroelectric- 
ity in this paper. 

We adopt the 1/4- filled extended Hubbard Hamilto- 
nian for holes on the two-dimensional anisotropic trian- 
gular lattice, which includes the degrees of freedom in a 
dimer. It is given by 

H = {Pn,mPn,m + V n ^ m n n n m ) + J Vt r V^( 1 ) 

<n.m> n 

where < n,m > denotes the neighbor site pairs, c„ i(T 
(cjj CT ) is the annihilation (creation) operator for a hole of 
spin a at the site n, fi n ,m {Vn,m) is the transfer integral 



FIG. 1: (Color online) Anisotropic triangular lattices for k- 
(BEDT-TTF) 2 X. 



(the Coulomb interaction energy) between the sites n and 



G n(J c m>a ), U is the on-site 



Coulomb interaction energy, n n , 



= r1 



and n n = 



Y^a n n,a- The 4x4 cluster of the system size N = 16 
with the periodic boundary condition is used. 

We calculate the quantum state excited by light 
pulse expressed by the vector potential A(t) = 
eAexp(—(t/D) 2 ) cos(w p i) of a Gaussian profile. We con- 
sider the one-photon excited states in this paper. To the 
first order in A, the solution \ip(t)) of the Schrodinger 
equation subject to the pulse is given by \ip(t)) = \ipo) + 
I^M 1 )^)). The ground state j^o) an d the first order solu- 
tion ^(-^(t)) are numerically calculated. 

We consider the first-order changes of the hole den- 
sity p { n\t) = ((^oKIV^H*)} +c.c.)/A and the bond 
order p^j m (t) = ((^ |p„, m |V (1) W) + cc.)/A induced by 
the pulse excitation. Since the absolute value of the bond 
order shows the strength of the bond, we define the dif- 
ference in the absolute value of the bond orders between 

\lp(t)) and \tpo) as \p\n}m{t) = \p n ,m + Pn,L(t)\ - \pn,m\, 

where p n . m is the bond order for the ground state. The 
bond is stronger in \tp(t)} than in \ipo) for IpI^L > 0. 

As shown in Fig. [H there are four nonequivalent bonds, 
and they are labeled by bl, b2, p, and q according to 
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Mori et al. 17 1. The bl bond is much stronger than the 



other ones, and the two BEDT-TTF molecules connected 
by the bl bond form a dimer. There are two dimers in 
the unit cell, and they are labeled by I and II as shown 
in this figure. The transfer integral /3bi and Vbi denote 
/3„ iin and V n , m for the bl bond, respectively, and the 
others are denoted in the same way. The transfer inte- 
grals for bl, p, and q bonds are fixed to the following 
values for K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br deduced from 
the extended Hiickel calculation [l7j ]: /?bi = 0.265 eV, 
(3 p = 0.109 eV, and /3 q = -0.038 eV. As for the Coulomb 
parameters, we use the following values: U = 0.7 eV 
El El, y bi = 0.45 eV, V h2 = V p = Vq = 0.25 eV. With 
these fixed parameters, we show the ^b2 dependence of 
various physical quantities in this paper because the spin 
frustration can be controlled by changing /3b2- 

We first show the physical properties of the ground 
state. The spin correlation function discontinuously 
changes at /3 b2 = 0.073 eV. For /3 h2 < 0.073 eV, the spin 
correlation is always negative between the dimers I and 
II, showing that the ground state is an antiferromagnetic 
(AFM) state. For /3 b2 > 0.074 eV, the spin correlations 
are much smaller than those for the AFM ground state, 
showing that the ground state is a paramagnetic (PM) 
state. Decreasing /?bi, the ground state discontinuously 
changes from the AFM state to the PM state. The PM 
ground state changes continuously as /?bi is decreased 
even to the value of /3b2 or /3 p . This shows that the PM 
state is stable without dimerization, and therefore it is 
not a dimer Mott insulator of a spin liquid state. We also 
calculate the probability W s that a dimer is singly occu- 
pied. The probability W s also discontinuously changes at 
the phase boundary, and W s > 0.85 holds all through the 
AF ground state region. This indicates that the AFM 
state is a dimer Mott insulator. On the contrary, W s 
decreases continuously to about 0.5 as /?bi is decreased. 
This result strongly suggests that the PM state is metal- 
lic. In the realistic parameter region, the ground state is 
either the AFM dimer Mott insulator or the PM metal, 
and the phase transition between them occurs by chang- 
ing various parameters. This result is consistent with 
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FIG. 2: (Color online) Light absorption spectrum a{ui) when 
the light is polarized (a) to the a-direction and (b) to the in- 
direction for various values of /3b2-The interval of /3b2 is 0.006 
eV. 
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FIG. 3: (Color online) The charge and bond order induced by 
the THz light polarized along the c-direction for /3b2 =0.068 
eV. (a) Nonequivalent bonds and sites are distinguished by 
colors and kinds (solid or dashed) of lines, (b) Pn (t) and (c) 



previous ones TH 19j]. The physical properties of these Hn,'™W are shown by the lines of the colors and kinds that 



AFM and PM states are basically unchanged in the re- 
alistic parameter range, and the essentially same results 
are obtained. 

In Fig. [21 we show the /3b2 dependence of light absorp- 
tion spectrum a(u) in the energy region of THz light. 
The dominant peaks in a(w) are in the region uj > 0.2 eV. 
In addition to these large peaks, there are much smaller 
peaks in the THz region both in the cases of light excita- 
tion from the AFM and PM ground states. It should be 
emphasized that such low energy peaks are not observed 
in the simple Mott and charge ordered insulators. 

We next show the time dependence of physical quan- 
tities induced by a THz light pulse. In all the following 
cases, the center frequency w p of the pulse is chosen to be 



specify the site n and the bond n, m, respectively, in (a). 



resonant to the largest peak in the energy region shown 
in Fig. [2j The used duration D = 100 eV -1 is large 
enough that the contributions from the large peaks in 
the higher energy region are negligible. 

We first show the results in the case of the excitation 
from the AFM ground state. In Fig. 02 we show the 
time dependence of charge and bond order induced by 
the THz light when the light is polarized to the c-axis at 
Pb2 = 0.068 eV. The bl, b2, p and q bonds are shown 
by the red, the green, the dark blue, and the light blue 
lines, respectively, in Fig. [3] (a). Moreover, there are 
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FIG. 4: (Color online) The schematic representation of bond 
order alternation between dimers. The ellipsoids represent 
the dimers, and the solid (dashed) lines show the stronger 
(weaker) bonds connecting the dimers. The solid (dashed) 
circles represent the charge-rich (charge-poor) sites. 



two nonequivalent p (q) bonds with different bond orders 
in the photoexcited state \ip(t)). These two nonequiva- 
lent bonds are distinguished by the solid and the dashed 
lines. Two sites forming a dimer also become nonequiva- 
lent, and they are distinguished by the solid red and the 
dashed green ellipsoids. We show Pn\t) and \p\n}m{t) in 
Figs. [3J (b) and (c), respectively, for the sites and the 
bonds shown in Fig. [3J (a). As seen from these figures, 
Pn\t) and IpI^L^) oscillate coherently. The quantities 
shown by the solid lines and the dashed lines oscillate out 
of phase. 

Before analyzing the present results, we introduce 
bond order alternation between dimers in much simpler 
ID quarter-filled electron-lattice coupled systems. The 
dimer Mott insulator with alternating bond lengths and 
bond orders between the dimers is schematically shown in 
Fig. The site connected by the stronger bond becomes 
charge-rich because electrons are more stabilized on the 
stronger bond. The state becomes the ground state with 
the deformed lattice in a parameter range US-HI. 

The collective motion induced by the THz light pulse 
can be regarded as the bond order alternation between 
the dimers. As seen from Fig. [3J (c), the amplitudes 
of the induced bond order oscillation for the p and q 
bonds are much larger than those for the other bonds, 
and therefore we focus on the p and q bonds. As seen 
from Fig. [3J (a), the interdimer bond orders alternate 
along the a-direction, and the amplitude of the alter- 
nation oscillates with time. Since charges are attracted 
by stronger bonds, when the bonds shown by the solid 
(dashed) lines become stronger, the sites shown by the 
solid red (dashed green) ellipsoids become charge-rich as 
seen from Fig. [3J (b) . In this way, the bond order and 
the charge density oscillations are strongly coupled in the 
collective mode, and the energy of the charge oscillation 
is significantly reduced by the coupling. The excitation 
induces the polarization along the c-direction due to the 
electric dipole moment on each dimer. 

It should be emphasized that the lattice deformation 
is not considered here. In the previous works, interdimer 
bond order alternation occurs as a result of electron- 
lattice coupling. On the other hand, in the present case, 
the induced collective motion of the bond order alter- 
nation coupled with polarization has purely electronic 
origin. We have found this purely electronic low-energy 
excitation also in the one-dimensional dimer Mott insu- 
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FIG. 5: (Color online) The charge and bond order dynamics 
induced by the THz light polarized along the a-direction for 
/3 b 2 =0.098 eV. (a), (b), and (c) The same as Fig. [3] 



lators without electron-lattice coupling. 

When the light is polarized to the a-axis, the collec- 
tive motion of the interdimer bond order alternation is 
induced by the THz light pulse. The collective motion is 
just like the case of the excitation from the PM ground 
state schematically shown in Fig. [5] (a). The bond or- 
ders for the p and q bonds alternate along the c-direction. 
The bond orders and the electric dipole moments oscil- 
late coherently in the mode, and the electric polarization 
along the a-direction is induced by the excitation. 

We next show the results in the case of the excitation 
from the PM ground state. When the light is polarized 
to the c-axis, the collective motion shown in Fig. [3J (a) is 
induced just like the case of the AFM ground state. 

In Fig. [5j we show the result when the light is polar- 
ized to the a-axis at /3b2 = 0.098 eV. As seen from this 
figure, the amplitude of the induced oscillation for the b2 
bonds is the largest between the bonds in this case, and 
it is much larger than the largest ones in the other cases. 
However, the induced dipole moment is much smaller in 
this case compared with the other cases. This can be un- 
derstood as follows. As seen from Fig. [5] (a), a dimer is 
connected by two solid b2 bonds or two dashed b2 bonds. 
The dominant motion of the b2 bonds does not accom- 
pany the bond order alternation between the dimers, and 
the coupling with the charge motion is very weak. 

As a result of electronic degrees of freedom in a dimer, 
weak charge disproportion is coupled with the interdimer 
bond order alternation. Therefore, there exist no such 
low-energy charge excitations in the Mott insulators. 
Furthermore, these collective modes occur irrespective of 
the AF spin order and the Mott insulating phase. Spin 
frustration also does not play an important role for these 
modes. The formation of the dimer and strong correla- 
tion effects are essential for these modes. 

The ferroelectric phase is not obtained within the 
present model. However, the ground state becomes fer- 
roelectric by introducing electron-lattice coupling. For 
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example, if the dimers denoted by I (II) in Fig. Q] 
shift to left (right), then the transfer integrals for the 
p (q) bonds become /3 p + sA (/3 q — srA), where A > 
shows the magnitude of bond length alternation, s = 1 
(—1) for the bonds shown by the solid (dashed) lines, 
and r > accounts for the difference in the electron- 
lattice coupling. Since j3 p > and /3 q < 0, the sign 
for the coupling term is opposite. We consider the state 
\ip) = \/l — Cf\ipo)+Ci\tpi), where C\ is a real constant, 
|V>i) = (V' (1) (' r )IV' (1) (' r )) _1/2 |V' (1) (' r )) wh en the light is 
polarized to the c-direction, and jpl^L^) for the bonds 
shown by the solid lines become maximum at t = r. The 
energy eigenstate is responsible to the absorption 
peak in the THz frequency region. The bond orders for 
\tp) for the x bonds are p x (l — C 2 ) + sp' x Ci ^\ — C\ + 
(l/2)24'Ci, where x=p or q, p x = (ipo\p n ,m\to) , p' x = 
(•>Po\Pn,m\ipi) + c.c. and p" = 2(ipi\p n ^ m \ip 1 ). The energy 
expectation value for with the distorted lattice minus 
the ground state energy Eq without lattice distortion is 
given by AE(C X , A) = aC\ + PC\A + 7A 2 , to the sec- 
ond order in C\ and A. Here a = E\ — Eq, E\ is the 
energy eigenvalue for P = N p (p' p — rp' q ), N p is the 
number of p or q bonds, 7 = N p (l/2)(K p + K q r 2 ), and 
Kx is the renormalized spring constant for x bond. AE 
is minimized when C\ = — /3/(2a)A and then AE = 
(7 - /3 2 /(4a))A 2 . Therefore, when a < 4/3 2 /7, lattice 
distortion occurs and is hybridized to | ^0} • As a re- 
sult, the ground state have a finite electric polarization. 
In the present case, \(3\ is very large because the contribu- 
tions from all the p and q bonds are negative as a result 
of the collective nature of the excitation. It is estimated 
that K K ~ 500 eV 1 for a-(BEDT-TTF) 2 X @. Using 
this value and setting r = 1, we obtain 4/3 2 /7 = 0.03 
eV. Since a is close to zero near the phase boundary, the 
ferroelectric phase is realized near the phase boundary, 
which is confirmed by a numerical calculation. 

With using the effective Hamiltonian where the in- 
terdimer part is treated perturbatively, the ferroelectric 
ground state is obtained without lattice deformation [24| , 
and this is not consistent with our result. We consider 
that this is because the interdimer part cannot be treated 
perturbatively. The effective Hamiltonian for the strong 
correlation case is also proposed (2S| . 

Based on the present results, dielectric anomaly ob- 
served in k-(BEDT-TTF) 2 Cu 2 (CN) 3 may be explained 
in the following way. As the temperature T is decreased, 
the system becomes close to the phase boundary. As a re- 
sult, the domains of the ferroelectric phase are generated 
for T < 60K, which results in the broad peak with relaxor 
like relaxation in dielectric constant. Around T = 6K, 
the ground state becomes ferroelectric, and this is the 
origin of the anomaly of various physical quantities 0- 
Hzl • I n particular, the present result is consistent 
with the unexpected experimental result that anomalous 
lattice response is observed only along the a-direction 
[26l 27 1 . Similar scenarios are proposed also in Refs. 
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Motional spectrum narrowing of an intermolecular 
phonon mode is observed in k-(BEDT-TTF) 2 Cu 2 (CN) 3 
around T = 30K [28J, and this can be attributed to the 
collective modes presented in this paper. These modes 
are strongly coupled with the phonon mode of interdimer 
bond length alternation. This is consistent with the fact 
that the narrowing occurs only in one specific phonon 
mode. Narrowing does not occur for T < 5K. This 
is because the ground state becomes ferroelectric and 
these modes disappear in the temperature region. Fur- 
thermore, the correlation time of the fluctuation mode 
roughly agrees with the energies of these collective modes 
apart from the phase boundary. 

This research was supported by CREST, JST. 
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